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ABSTRACT 

We present an analysis of the stellar kinematics of the Large Magellanic Cloud based on ^5900 new 
and existing velocities of massive red supergiants, oxygen-rich and carbon-rich AGB stars, and other 
giants. After correcting the line-of-sight velocities for the LMC's space motion and accounting for 
asymmetric drift in the AGB population, we derive a rotation curve that is consistent with all of the 
tracers used, as well as that of published H I data. The amplitude of the rotation curve is vq = 87 ± 5 
km s _1 beyond a radius Ro = 2.4 ± 0.1 kpc, and has a position angle of the kinematic line of nodes 
of 9 = 142° ± 5°. By examining the outliers from our fits, we identify a population of 376 stars, or 
>5% of our sample, that have line-of-sight velocities that apparently oppose the sense of rotation of 
the LMC disk. We find that these kinematically distinct stars are either counter-rotating in a plane 
closely aligned with the LMC disk, or rotating in the same sense as the LMC disk, but in a plane that 
is inclined by 54° ± 2° to the LMC. Their kinematics clearly link them to two known H I arms, which 
have previously been interpreted as being pulled out from the LMC. We measure metallicities from 
the Ca triplet lines of ^1000 LMC field stars and 30 stars in the kinematically distinct population. 
For the LMC field, we find a median [Fe/H]=— 0.56 ± 0.02 with dispersion of 0.5 dex, while for the 
kinematically distinct stars the median [Fe/H] is —1.25 ± 0.13 with a dispersion of 0.7 dex. The 
metallicity differences provide strong evidence that the kinematically distinct population originated 
in the SMC. This interpretation has the consequence that the H I arms kinematically associated with 
the stars are likely falling into the LMC, instead of being pulled out. 

Subject headings: Magellanic Clouds - galaxies: kinematics and dynamics - galaxies: interactions 



1. INTRODUCTION 

The Magellanic Clouds have always been full of sur- 
prises, many of which stem from the fact that the Clouds 
are in the midst of interacting with each other and with 
the Milky Way. A famous surprise was the discovery 
of the Magellanic Stream trailing the Clouds (Wannier 
& Wrixon 1972; Mathewson et al. 1974), which is now 
known to have a companion Leading Arm (Putman et 
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al. 1998), and which taken together extend for ^200° 
across the sky (Nidever et al. 2010). Other striking ev- 
idence of the interaction is found in the Small Magel- 
lanic Cloud's and the Large Magellanic Cloud's intrin- 
sically elongated shapes (Caldwell & Coulson 1986, van 
der Marel 2001), the possible warps of the LMC (Olsen & 
Salyk 2002, Nikolaev et al. 2004, Subramaniam & Subra- 
manian 2009), and the young population of stars formed 
from tidal debris in the Magellanic Bridge (Harris 2007). 
A recent surprise is that with new proper motion mea- 
surements of the Clouds (Kallivayalil et al. 2006a, 2006b; 
Piatck et al. 2008), the trajectories of the Clouds indi- 
cate that they are likely on first or second approach to 
the Milky Way (Besla et al. 2007, Boylan-Kolchin et al. 
2010), such that most of the interaction effects, e.g. the 
Stream, must be due to the LMC-SMC binary pair in- 
teraction, not an interaction with the Milky Way (Besla 
et al. 2010). Many further surprises are likely in store, as 
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foreshadowed by Munoz et al. (2006), who found LMC 
stars 20 kpc away in the foreground of the Carina dwarf, 
and Saha et al. (2010), who showed that the LMC main- 
tains an exponential disk profile out to at least 12 disk 
scale lengths. 

An observed feature found in some merged or inter- 
acting disk galaxies is the presence of extended counter- 
rotating stellar disks. The first, remarkable discovery 
of such a system is NGC 4550, an SO galaxy with two 
roughly equal mass counter-rotating stellar disks (Rubin 
et al. 1992, Rix et al. 1992). The NGC 5719/13 system 
appears to be an example of an interacting system where 
a significant counter-rotating stellar disk is forming (Ver- 
gani et al. 2007). However, Kuijken et al. (1996) found 
no additional counter-rotating cases in a sample of 28 
SOs, and concluded that no more than ^5% of the stars 
in these galaxies could be counter-rotating, a testament 
to the difficulty of detecting such stars. 

We report on the first results from an extensive spec- 
troscopic study of the LMC. As described below, we have 
found that the LMC contains a significant (~ 5%) stel- 
lar population whose kinematics indicate either counter- 
rotation in a plane closely aligned with the LMC disk, 
or rotation in a plane that is highly inclined to the 
LMC disk. Evidence for counter-rotation in the core 
(radius<2.5 kpc) of the LMC was previously presented 
by Subramaniam & Prabhu (2005), while Kunkel et al. 
(1997) presented evidence for an out-of-plane polar ring 
in the LMC. Here, we present clear evidence of such a 
kinematically distinct stellar population, as well as evi- 
dence that the population is kinematically cold and orig- 
inated in the SMC. The LMC thus appears to have ac- 
creted stars from the SMC, and may be a unique labora- 
tory for studying the formation counter-rotating disks. 

2. OBSERVATIONS 

The primary observational data used in this paper was 
collected during the nights 20 - 26 Nov 2007 at the CTIO 
4-m Blanco t elescope with the Hydra- CTIO multi-fiber 
spectrograph (jBarden fc Ingersonll 19981 ) , with some addi- 
tional data obtained on 19 Nov 2007 during time shared 
with instrument checkout and engineering. All of the 
nights were clear, with only some thin cirrus on the 
night of 21 Nov. We used Hydra-CTIO with the 400mm 
Schmidt camera and SITe 2Kx4K CCD, giving a gain 
of 0.84 e~/ADU and 3 e~ read noise. We observed with 
the red-blazed KPGLF grating at a central wavelength of 
8000A and an OG515 order-blocking filter. We masked 
the fibers with a 100/iin slit mask, yielding a 2-pixel 
spectral resolution of R ~ 7000 and velocity precision of 
~2— 3 km s , over a wavelength range of ~6850-9150A. 
In addition to our program targets, each night we ob- 
served several spectral template stars from a list roughly 
matching the spectral types of our targets; these tem- 
plate objects consisted of the stars HD223311 (spectral 
type K4III), HD16115 (R-type carbon star), HD80170 
(K5 III), and BN Mon (N-type carbon star). 

Our targets for the Hydra-CTIO observations were 
selected using the Spitzer SAGE survey of the LMC 
(M eixn er et al.l 120061 ) in combin ation with the near- 
infrared photometry from 2MASS (jSkrutskie et al.ll2006D 
and optical photome try from the Magel lanic Clouds 
Photometric Survey (|Zaritskv et al.l 120041) . We iden- 



tified candidate LMC red supergiants (RSGs; evolved 
from stars with initial masses >10 M©), giants evolved 
from stars with initial masses ~5 M Q , carbon-rich AGB 
stars, oxygen-rich AGB stars, as well as "extreme" 
AGB stars. The definitions of these classes were taken 
from the analysis of the Spitzer SAGE photometry by 
iBlum et "all ()2006f ) and by matching the above photo- 
metric databa ses with the pos i tions of previously con- 
firmed RSGs (iMassev k. Olsenl l200"3l ) and carbon stars 
(|Hardv et alJ 120011 ). We chose these target classes both 
because they span a range of stellar ages and because 
their numerous spectral features provide precise veloc- 
ity measurements. We also chose to observe 558 candi- 
date bright LMC main sequence stars; nearly all of these 
turned out to lie in the Milky Way foreground, however, 
and are not discussed further. We assigned the chosen 
targets to individual Hydra-CTIO fibers using the pro- 
gram hydraassign (Massey 1995), which resulted in 62 
configurations with ~95 targets per configuration. All of 
the remaining ~30 fibers were used to measure the sky. 
Fig. 1 shows the J — [3.6], [3.6] color-magnitude diagram 
(CMD) of stars for which we measured spectra, com- 
pared to all of the available sources. In all, we obtained 
useful spectra of 4734 targets, which excludes the 558 
main sequence stars in the Milky Way foreground that 
we observed. 4567, or 96%, of these targets were found 
to have velocities consistent with LMC membership. 

3. DATA REDUCTION AND ANALYSIS 

We processed our data using IRAB3 routines. In brief, 
we used CCDPROC to subtract the overscan region, to 
subtract a bias frame, and to trim extraneous pixels from 
the image edges. W e then used the routine L. A. Cosmic 
(|van Dokkuml l2001) to remove cosmic rays from the in- 
dividual spectral images. We finally used DOHYDRA to 
extract the individual fiber spectra, apply the dome flat 
field, solve for and apply the dispersion solution based on 
our He/Ne/Ar arc lamp observations, and subtract the 
average sky spectrum. 

We measu red velocities for the L MC stars using cross- 
correlation (jTonrv &; Davisl Il979f ) with our template 
spectra of stars with known radial velocities, using the 
IRAF task FXCOR. We used FXCOR to first apply a 
ramp filter in the Fourier domain to dampen the highest 
frequency (e.g. noise) and lowest frequency (e.g. con- 
tinuum) Fourier components, and then to fit the narrow 
cross-correlation peaks with Lorentzians to find their cen- 
troids. We performed the cross-correlation in wavelength 
windows that avoid the strongest telluric absorption fea- 
tures, and used the template with the strongest corre- 
lation peak, as measured by the r parameter, to make 
the velocity measurement. We used the errors reported 
by FXCOR, which derive from r, as the velocity error. 
We accepted only those measurements with r > 5, which 
in practice translated to velocity errors <5 km s _1 ; the 
average velocity error was measured to be 2.2 km s _1 . 
We confirmed the stability of our velocity measurements 
by cross-correlating template spectra taken at different 
times during each night and on different nights. In all 

3 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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cases, the template velocities agreed to within 1—2 km 
s _1 , such that velocity drifts across nights should be no 
larger than the measured random errors of —2 km s _1 . 

A possible concern with regard to the velocities is sys- 
tematic error introduced by the mismatch in the spec- 
tral types of our program stars with that of our template 
stars. In particular, a number of our program stars have 
M-type spectra, whereas our closest velocity templates 
for these stars have mid-K-type spectra. Despite the vis- 
ible differences between K giant and M supergiant spec- 
tral types (most notably the strong TiO bands present 
in the M supergiants), there are still many spectral lines 
common to both types, e.g. the Ca II triplet lines at 
-8600 A. In Fig. 2, we show representative spectra of a 
carbon star, an oxygen-rich AGB star, a red supergiant, 
and a ~5 M© giant compared to the velocity template 
used for each star. The spectra are not flux-calibrated. 
The color of the carbon star suggests that it is R-type, 
whereas the template star shown, BN Mon, is N-type. 
Judging from the prominent TiO absorption at ~7100 
A and ~8400 A, the O-rich AGB star is a mid-M giant 
(Silva & Cornell 1992), while the template shown is type 
K4III. The presence of TiO absorption in the spectrum 
of the red supergiant suggests that is a late K supergiant, 
whereas the template is K5III. Finally, the ~5 M© giant 
appears to be a mid-K-type star, as is its template. As 
seen in the figure, we obtained strong, symmetric cross- 
correlation peaks for each of the spectra, even in cases 
where our targets had later spectral types than the tem- 
plates. In addition, after applying the computed velocity 
shift, the 8662 Ca II line clearly falls at the same wave- 
length in the template and the program star. We thus 
do not think that template mismatches are a significant 
source of systematic error in the velocities. 

As a final check, we compared our velocity measure- 
ments of stars that overlapped with previous work. For 
the 19 RSGs in common with Massey & Olsen (2003), we 
found an average difference of Uhci ~~ u hoi,MO=3.9 ± 1.2 
km s , where the error is the standard error of the 
mean, and a median difference of 2.8 km s _1 . For the 97 
RSGs in common with Prevot et al. (1985), who observed 
their sample with the precise velocity scanner CORAVEL 
(Baranne et al. 1979), we found an average difference of 
Vhci — fhci,Prcvot=0.6±0.5 km s _1 . For the 91 carbon stars 
in common with Hardy et al. (2001), the average differ- 
ence is Vhei — «hei,Hardy =— 2.4±0.4km s" 1 . The very close 
agreement with previous work reinforces our conclusion 
that systematic errors in our velocities are small. 

We performed our kinematic analysis on the 4567 LMC 
RSGs, giants, and AGB stars observed with Hydra- 
CTIO, the 857 carbon stars observed by iKunkel et al.l 
((1991 and Hardy et al. (2 01), a nd t he 481 RSGs 
observ ed by iMassev fc Olsenl ((2003) and iPrevot et al.1 
(1985)-a total of ~5900 stars. As noted in the com- 
prehensive paper of van der Marel et al. (2002; hereafter 
vdM02), interpretation of radial velocities in the LMC 
must account both for contribution of the LMC's space 
motion projected into the line of sight, as a function of 
position on the sky, and for the possibility of prcces- 
sio n and nutation of the LM C disk. Indeed, as shown 
bv Ivan der Marel et al.l (|2002f l. the space motion contri- 
bution has a larger effect on the line-of-sight velocity 
distribution than do the LMC's internal motions. The 



accurat e LMC prope r measurement bv (K allivavali l et al.l 
(|2006al) . updated bv lPiatek etaLl (|2008l ). fixes the space 
motion contribution and thus greatly improves the accu- 
racy of the kinematic analysis. 

As we did in Olsen & Massey (2007, hereafter Paper 
I) , we followed the procedure outlined by vdM02 to solve 
for the LMC's systemic velocity v sys , the rate of change 
of the disk's inclination di/dt, the position angle of the 
kinematic line of nodes 9, and the LMC's internal ro- 
tation curve v(R). We fixed the location of the LMC's 
dynamical center to the position determined by vdM02, 
a'™ = 5 h 27 m 36 s , 6 C m = -69°52'12", as this position 
coincides closely with the LMC's geometric center as de- 
termined from structure out to large radii (van der Marel 
& Cioni 2001). We furthermore set the inclination of the 
disk to i — 34?7, as determined by studies of the LMC's 
geometry (van der Marel & Cioni 2001), and parameter- 
ized the LMC's rotation curve by a linear function with 
a value of km s" 1 at radius R = and rising to a value 
vq at radius Rn , with v(R) — vq for R > Rq (as done 
bv lPiatek etafl ((20081 )). We calculated errors in the fit- 
ted quantities by creating 10000 Monte Carlo samples of 
our data, refitting the parameters, and calculating the 
dispersion in the parameters. The Monte Carlo samples 
were produced through a bootstrap process, in which we 
populated the samples by randomly drawing stars from 
our original dataset, but allowing any individual star to 
be drawn any number of times in the same sample. 

Fig. 3 (left) shows a representation of the result of our 
fit to the 738 RSG line-of-sight velocities, which yielded 
parameters of v sys = 263 ± 2 km s _1 , di/dt—— 184 ± 81° 
Gyr-\ 9 = 142 ± 5°, v = 87 ± 5 km s~\ and R„ = 
2.4 ± 0.1 kpc. Following Paper I, in making this figure 
we first removed the contribution of the LMC's space 
motion and the contribution from di/dt, then assumed 
that the remaining velocity signatures were due to in- 
plane circular motions in the LMC disk. For comparison, 
we a lso show th e emission- weighted mean H I velocitie s 
from lKim etai] ((19981 ) and iStavelev-Smith et all (|2003l) , 
with the same corrections applied as to the RSGs. Note 
that while the fit was performed on the line-of-sight ve- 
locities, i.e. the projected velocities due to the LMC's 
internal kinematics and space motion, in this figure we 
show the de-projected circular velocities implied by our 
fit. For stars close to the line of nodes, the factor needed 
to deproject the velocities becomes large, such that we 
have excluded stars with projection factors larger than 
5 from the left panel of Fig. 3 to avoid confusion in the 
visual presentation; the H I gas in the region containing 
these excluded stars is shown as grayscale in the right 
panel of Fig. 3. 

As in Paper I, Fig. 3 shows that much of the H I gas de- 
fines a flat rotation curve; in contrast to Paper I, however, 
the much larger sample presented here yielded a new so- 
lution that shows few significant RSG velocity outliers 
compared to the H I, and defines what appears to be a 
nearly identical rotation curve to the H I. We also iden- 
tify the velocity signatures of the H I arms S, E, and B 
(jStavelev-Smith et al.ll200l and E2 (Paper I), and their 
spatial location in the right-hand panel of Fig. 3. We 
confirm the Paper I result that few if any RSGs are as- 
sociated with the H I arms, with a larger sample of RSG 
velocities. 
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In Fig. 4 (left), we compare the implied in-plane cir- 
cular velocities of our entire sample, which is predom- 
inantly made up of carbon-rich and oxygen-rich AGB 
stars, with the kinematic model ht to the RSG-only sam- 
ple. The bulk of the stars trace the envelope of the flat 
rotation curve, but with many outliers with lower implied 
circular velocities, including many negative circular ve- 
locities. Examining the line-of-sight velocity residuals of 
the full sample about the fit to the RSGs (Fig. 4 right 
panel) reveals an offset in the peak of the distribution, 
< uios.fuii ~ uios.RSG >= -4.2 km s~\ and a dispersion 
of 26 km s _1 . The offset can be explained as asymmetric 
drift, with the magnitude consistent with that expected 
from the observed velocity dispersion, following vdM02. 
In Paper I, we found substantially different rotation ve- 
locities for LMC RSGs, H I, and carbon stars. With 
the much larger sample of stars studied here, it appears 
that these differences are only as large as expected by 
asymmetric drift, such that all populations have basic 
kinematics that are consistent with each other. 

To understand the outliers better, we selected all stars 
with Afi os > 50 km s -1 (which have velocities in excess 
of those predicted by our kinematic model) and plotted 
them as red points in the left panel of Fig. 4, as well as all 
stars with Avi os < —50 km s" 1 , which are plotted as blue 
points. Strikingly, the majority of these outliers fall in 
the region where we see the kinematic signatures of the E 
H I arm, and its E2 extension, and the B arm. In Paper I, 
we observed carbon stars in these regions of the plot, and 
interpreted them as stars that were being tidally stripped 
from the LMC along with the H I gas; indeed, many of 
the stars appeared co-located with the H I arms. With 
our larger dataset, we see many more of these stars (376), 
as well as roughly equal numbers with positive (~40%) 
and negative (~60%) velocity residuals. Fig. 5 shows the 
location of the stars with |Ai>i os | > 50 km s _1 on the sky, 
where the colors indicate the line-of-sight velocities after 
correction for the LMC's space motion and for di/dt. 
The underlying rotation of the LMC is visible in the bulk 
of the stellar population; it is worth noting that the near 
side of the LMC lies in the northeast (Caldwell & Coulson 
1986), such that the sense of rotation is clockwise on the 
page. The stars with |Awi os | > 50 km s _1 , however, not 
only have large residuals from the model, but clearly have 
line-of-sight velocities that oppose the sense of rotation of 
the disk. Almost all of them also lie inside the outermost 
radius at which the LMC rotation curve can be traced. If 
these stars lie in the LMC plane, then they must indeed 
be counter- rotating compared to the LMC. If they lie 
in a structure with a very different geometry, e.g. in a 
plane whose near side is in the southwest rather than the 
northeast, then their sense of rotation could be the same 
as the LMC disk. In either case, we think it is untenable 
to maintain the interpretation that these kinematically 
distinct stars have been tidally stripped from the LMC, 
as we concluded in Paper I based on a smaller sample 
of stars and on our expectation that the H I arms with 
which the stars were associated were being pulled out of 
the LMC. We note that the fraction of this population, 
>5%, is similar to the 7% of carbon stars that Graff 
et al. (2000) found comprised a kinematically distinct 
population in their analysis. The similar fractions may 
be fortuitous, however, as we see no clear sign of the 



kinematically distinct population in the raw histogram 
of heliocentric velocities, contrary to what was seen by 
Graff et al. with their smaller sample. 

To examine the geometry of the kinematically distinct 
population, we made the assumption that the stars are on 
circular planar orbits in the LMC's potential, for which 
we know the rotation curve from the RSGs above. Fix- 
ing i> sys , vq , and Rq at the values determined for the 
LMC RSGs, we fit for the rate of change of inclina- 
tion (di/d£)KDP, the angle of the kinematic line of nodes 
Okdp, and inclination ikdp of the kinematically distinct 
population. We found two solutions, the first being a 
counter- rotating solution with (di/di)KDP = 290° ± 50°, 
&kbp = 177 ± 7°, and *kdp = 20 ± 3°. This solu- 
tion would indicate that the kinematically distinct pop- 
ulation lies in a plane that is twisted by ~ 35° com- 
pared to the LMC disk and with an inclination that is 
~ 15° shallower with respect to the line of sight, but is 
quantitatively speaking roughly coplanar with the LMC. 
This solution is similar to that adopted by Subrama- 
niam & Prabhu (2005), who modeled the LMC core as 
two disks with identical inclination but lines of nodes 
twisted by 40° from each other. The second solution 
has the stars rotating in a clockwise sense on the sky, 
but with (di/dt) KDP = 300° ± 50°, 6» KD p = 175 ± 7°, 
and Ikdp = — 19 ± 2°. In this solution, the kinemat- 
ically distinct stars lie in a plane that is inclined by 
~ 55° compared to the LMC disk, similar to the frag- 
mented polar ring model of Kunkel et al. (1997). Nei- 
ther solution depends strongly on the fitted value of 
(dz/dt)xDp; fixing (di/di)i<DP equal to the value derived 
for the LMC RSGs led to solutions with nearly identi- 
cal zkdp, but with a more highly twisted line of nodes, 
#kdp = 217°±7°. Both solutions have line-of-sight veloc- 
ity dispersions about the solutions of 15 km s _1 , which is 
significantly smaller than the raw observed velocity dis- 
persion of 26 km s _1 for the kinematically distinct popu- 
lation. The small velocity dispersions about the solutions 
implies vja ^5 for the kinematically distinct population, 
demonstrating that the population is kinematically cold. 

We next gathered evidence to understand the origin of 
the kinematically distinct population. Fig. 6a shows the 
J — [8.0], [3.6] CMD for the kinematically distinct popula- 
tion compared to all of the available sources. The figure 
shows that nearly all of the kinematically distinct stars 
are O-rich and C-rich AGB stars, but are otherwise not 
remarkably different from the normal LMC AGB pop- 
ulation, suggesting that they are at the same distance 
as the LMC. The only difference appears to be that the 
normal LMC population contains some luminous O-rich 
and C-rich AGB stars that are not seen in the kinemati- 
cally distinct population. To quantify this difference, we 
counted the number of LMC AGB stars found within the 
red outline in Fig. 6a, which are brighter than those in 
the polygon that contains the majority of the kinemat- 
ically distinct population. If the kinematically distinct 
population were drawn from the normal LMC AGB pop- 
ulation, we would expect to observe 10±3 AGB stars 
brighter than those in the polygon, whereas we observe 
one star at the edge of the outline, with [J-8.0]^2.5 and 
[3.6]^7. This star could be a bona fide bright AGB star 
or a red LMC RSG with a peculiar velocity. A much more 
striking difference, shown in Fig. 6b, is seen in the spa- 
tial distribution of the kinematically distinct population 
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compared to the normal LMC AGB population. While 
normal LMC AGB stars are found at all locations within 
the SAGE survey, they are most heavily concentrated in 
the LMC Bar. The kinematically distinct population, 
however, clearly avoids the Bar, implying a different ori- 
gin from that of the normal LMC AGB stars. 

We used our spectra to measure metallicities from the 
Ca II near-infrared triplet lines for a subset of stars in 
the kinematically distinct population and in the normal 
LMC population. We used only those spectra without 
carbon star features and with S/N per resolution ele- 
ment in excess of ~20, which yielded measurements for 
994 normal LMC stars and 30 stars in the kinematically 
distinct population. We followed Cole et al. (2004) to fit 
summed Gaussian and Lorentzian profiles to the spec- 
tral lines, integrated these fits to derive their equivalent 
widths, and formed Ca indices by summing the equiva- 
lent widths of the triplet lines. After using the optical 
photometry and extinction map of Zaritsky et al. (2004) 
to correct the V magnitudes of our stars for extinction, 
we used the calibration of Cole et al. (2004) to derive 
metallicities from the Ca indices and V — Vhb, the mag- 
nitude differences between our stars and the LMC hori- 
zontal branch. Following Cole et al. (2005), we adopted 
Vhb = 19.22 for all of the stars. Note that although most 
of our stars are brighter than the globular cluster stars 
used to calibrate the Ca index/U — VnB/[Fe/H] relation- 
ship, we found that the normal LMC population defines 
the same slope in the Ca index/1/— Vhb diagram as found 
by Cole et al. (2004), suggesting that the same calibra- 
tion applies. For the 994 normal LMC stars, we found 
a metallicity distribution with a peak at [Fc/H]=— 0.45, 
a median of [Fe/H]=— 0.56 ± 0.02, and a dispersion of 
0.5 dex compared to a median error of 0.15 dex. These 
measurements are in good agreement with the LMC Bar 
metallicity distribution measured by Cole et al. (2005), 
who found a median [Fe/H]=— 0.4 with a dispersion of 
0.3 dex. For the 30 stars from the kinematically distinct 
population for which we could measure the Ca triplet 
lines, we instead found a median [Fe/H] of —1.25 ± 0.13 
with dispersion 0.7 dex compared to a median error of 
0.17 dex. 

4. DISCUSSION AND CONCLUSIONS 

Where did the kinematically distinct population orig- 
inate? The stars have a ki nematic signature that links 
them to H I arms E and B. IStavelev-Smith et al.l (2003) 
identify these arms as connecting to the Leading Arm 
(Putman et al. 1998) in the case of arm E and the Mag- 
ellanic Bridge and the SMC in the case of arm B. The 
most tempting and obvious answer is thus that the stars 
originated in the SMC, and have fallen into the LMC as 
a result of the LMC-SMC-Milky Way interaction. How- 
ever, both arms E and B are thought to be outflows, 
perhaps driven by a combination of star formation ac- 
tivity and tidal forces (Nidever et al. 2008), whereas an 
association with stars captured from the SMC would ap- 
pear to require that they represent gas infall. We cannot 
discount the possibility of infall, as without knowledge 
of the gas geometry, its direction of motion is inherently 
ambiguous. An infall hypothesis for the arm E and B gas 
would indeed yield a natural explanation for two puz- 
zling aspects of the LMC. First, H I maps show that the 
southeastern rim of the LMC, near the location of 30 Dor, 



contains an overdensity of H I gas (e.g. Kim et al. 1998), 
and a similar overabundance of CO and giant molecular 
clouds (Fukui et al. 1999). The morphology of the gas 
led de Boer et al. (1998) to conclude that the overabun- 
dance exists because gas piles up on the leading edge of 
the LMC due to ram pressure as the LMC moves through 
the Galactic halo, while Nidever et al. (2008) identify the 
overdensity as the source of the outflows which they con- 
cluded feed arms E and B. In the infall hypothesis, the 
overdensity would instead be produced as the infalling 
gas encounters the LMC disk, either in a nearly head- 
on collision (if the kinematically distinct population is 
counter-rotating) or in a "T-bone" collision (if the ge- 
ometry of the kinematically distinct population is highly 
inclined with respect to the LMC). The stars would con- 
tinue through the site of the collision unimpeded, while 
the gas would shock as it encountered the LMC gaseous 
disk. The infall hypothesis could thus account for a sec- 
ond puzzling aspect of the LMC, namely the unique prop- 
erties of the 30 Dor star forming complex. As discovered 
by Dolphin & Hunter (1998), the region of Constellation 
III in the LMC formed from roughly the same gas mass 
as 30 Dor, but 30 Dor has yielded several times more stel- 
lar mass and at higher spatial concentration, including 
the dense central cluster like R136. Dolphin & Hunter 
suggested that this was because the initial distribution of 
gas in 30 Dor was much more compact, something which 
a violent collision of gas would naturally explain. 

Our measurements of the Ca triplet abundances pro- 
vide a check on the hypothesis that the kinematically dis- 
tinct stars came from the SMC. De Propris et al. (2010) 
used the Ca triplet to measure abundances of hundreds of 
SMC red giants at the periphery of the SMC, and found 
a metallicity distribution with mean [Fe/H]=— 1.35±0.10 
and a dispersion of 0.65±0.08 dex, in very good agree- 
ment with our measurement of the metallicity distribu- 
tion of the kinematically distinct population. Although 
we would like to confirm this result with high signal- 
to-noise abundance measurements of many more stars, 
the Ca triplet abundances provide strong evidence that 
the kinematically distinct population did originate in the 
SMC. We find further evidence that the stars originated 
in the SMC from Fig. 7, which shows a comparison of the 
J - [3.6], [3.6] CMD of the SMC (Gordon et al., submit- 
ted) with the CMD of the kinematically distinct stars as 
shifted to the distance of the SMC. In making this fig- 
ure, we have assumed that the SMC has distance mod- 
ulus 0.4 mags larger than the LMC, which is based on 
an LMC distance modulus of 18.5 and an SMC modulus 
of 18.9 (Storm et al. 2004). We have ignored extinction 
by dust, since for our stars we found < Ay >= 0.4, 
which translates to negligibly small Aj ~ 0.1 (Schlegel 
et al. 1998) and A [3 . 6] = 0.02 (Indebetouw et al. 2005). 
The figure demonstrates that the CMDs of the SMC and 
kinematically distinct LMC stars are qualitatively simi- 
lar, particularly if we restrict the SMC CMD to include 
only stars beyond radii of ~2 degrees from the optical 
center of the SMC, which we take to be a = h 49 m 47 s , 
5 = -72°53'40" (Westerlund 1997). At these outer radii, 
the most luminous AGB stars that are seen in the full 
SMC CMD, but not in the kinematically distinct popu- 
lation, no longer appear. To quantify this distinction, we 
counted the number of SMC AGB stars found within the 
area defined by the red outline in Fig. 6a (shifted verti- 



6 



cally by +0.4 magnitudes to account for the difference in 
distance modulus between the SMC and LMC) and com- 
pared it to the number of stars with properties similar to 
the kinematically distinct population, as defined by the 
polygon in Fig. 6a (again shifted for distance modulus). 
If the kinematically distinct population were drawn from 
the full SMC AGB population, we found that we would 
expect to observe 6±2 bright AGB stars in the sample, 
whereas we observe one. If the kinematically distinct 
population were instead drawn from the SMC AGB pop- 
ulation beyond a 2° radius from the center, we would 
expect to see 1±1 bright AGB star, and we observe one. 
Although these values are subject to small number statis- 
tics, such agreement with the stellar populations in the 
outer regions of the SMC is just what we would expect 
if the kinematically distinct stars were stripped from the 
SMC in a close passage between the Clouds. 

Coupled with our comparison of the stellar and H I 
kinematics, our results require a reinterpretation of the 
formation of H I arms E and B, and suggest that we 
should also examine our findings in light of models de- 
signed to form the Leading Arm feature and the Mag- 
ellanic Stream. A recent model is that of Besla et al. 
(2010), who have successfully reproduced the basic prop- 
erties of the Stream and the Leading Arm. In their 
model, the Clouds are a binary system that experienced 
a close interaction 1.2 Gyr ago, during which the Stream 
was formed from gas tidally removed from the SMC by 
the LMC. The SMC is modeled as having a compact 
stellar disk and an extended gas disk, resulting in no 
stripped SMC stars, and the Stream without an observ- 
able stellar component. Our result argues, however, that 
~ 5% of the stars in the LMC came from the SMC, such 
that the Besla et al. model needs some modification. In- 
deed, with a more extended SMC stellar disk, the Besla 
et al. model predicts that SMC stars should be captured 
by the LMC (Besla, personal communication); investi- 
gating the geometry and kinematic properties of these 
captured stars requires more modeling work. The exis- 
tence of accreted SMC stars in the LMC would then also 
predict that SMC stars should be present in the Stream. 
Their surface brightness would still be predicted to be 
lower than current detection limits, but should be reach- 
able by deep wide-field imaging with current or future 
facilities. 

In summary, we have used ~4600 new spectra and 
^1300 published velocities of primarily red supergiants 
and AGB stars to study the internal kinematics of the 
LMC. Our main conclusions are: 

1. After correcting the line-of-sight velocities for the 
LMC's space motion and accounting for asymmet- 
ric drift in the AGB population, we have found a 
rotation curve that is consistent with the red su- 
pergiants, the H I, and the AGB stars. 



AGB population. Moreover, they have a spatial 
distribution that avoids the LMC Bar, contrary to 
the normal LMC AGB population. 

4. Our simple calculation finds that the kinematically 
distinct stars are either counter-rotating in a plane 
closely aligned with the LMC disk, or rotating in 
the same sense as the LMC disk, but in a plane that 
is inclined by 54° ±2° to the LMC. Their kinematics 
clearly link them to two known H I arms, which 
have previously been interpreted as being pulled 
out from the LMC. 

5. The metallicities of the kinematically distinct pop- 
ulation, measured from the Ca triplet lines, are low 
compared to the LMC but in good agreement with 
the metallicities of SMC giants. 

6. The J — [3.6], [3.6] color- magnitude diagram of the 
kinematically distinct population is a good visual 
match to that of the outer regions of the SMC when 
shifted to the same distance. 

7. Our conclusion that the kinematically distinct stars 
found in the LMC were accreted from the SMC 
implies that the associated H I arms, previously 
thought to be flowing away from the LMC, are 
likely to be falling in. This interpretation could 
help to explain the overdensity of gas seen on the 
LMC's southwestern rim, and identifies a possible 
source of fuel for star formation in 30 Doradus. 
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2. We have discovered a population of stars, repre- 
senting ~5% of our sample, that has distinct kine- 
matics. 

3. The kinematically distinct population is comprised 
almost entirely of AGB stars, but do not have the 
number of bright AGB stars that would be ex- 
pected if they were drawn from the normal LMC 
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Figure 1. The color-magnitude diagram of all sources toward the LMC measured by SAGE (black points; Meixner et al. 2006). Filled 
circles are those stars for which we obtained spectra. The colors of the filled circles denote our rough classification of the stars, based on 
the criteria established by Blum et al. (2006): yellow for red supergiants, cyan for giants, green for oxygen-rich AGB stars, red for carbon 
stars, and black for "extreme" AGB stars. Boxes identify areas within which we randomly selected a star to represent its class. 
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Figure 2a. Spectra of a representative LMC carbon star (a) and oxygen-rich AGB star (b) are shown, along with the velocity template used 
for each star. Left column: The LMC spectra are shown in black, with grey marking wavelengths contaminated by terrestrial absorption. 
The spectra have been shifted in wavelength according to the velocity computed by cross-correlation with the template. The template 
spectra are shown in red, are labelled with names, and have been shifted vertically to avoid confusion. Middle columniThe LMC spectra 
are overlaid with the templates in a narrow wavelength region centered on the 8662 A Ca II triplet line. Right column: The correlation 
strength computed by FXCOR is shown as a function of the relative velocity shift. The upper panels show the correlations over their full 
range, while the lower panels present a zoomed-in view of the peaks at the computed relative shift. The grey horizontal lines mark the zero 
levels. 




Figure 2b. Spectra of a representative LMC red supergiant (c) and ~5 Mq giant (d) are shown, along with the velocity template used 
for each star. See Fig. 2a for a description of the plots. 
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In-plane radius (kpc) x (Degrees) 

Figure 3. The rotation curve of the LMC as fit to red supergiants. Left: The line-of-sight velocities of red supergiants (circles) and of 
H I gas (grayscale; Kim et al. 1998) have been deprojected to show their in-plane circular velocities, as described in the text. The green 
line is our parameterized fit of the rotation curve to the RSGs. The colored polygons identify the distinct kinematic components discussed 
in paper I, most of which were identified by Staveley-Smith et al. (2003): the flat rotation curve (red), arm S (purple), arm E and its 
E2 extension (yellow and green), and arm B (blue). Right: The H I gas and red supergiants are shown as projected on the sky, with 
color-coding corresponding to the kinematic components identified at left. 




In-plane radius (kpc) Av ]Ils (km s ') 

Figure 4. The kinematics of all stars studied in this work. On the left, the red supergiant rotation curve (white line) is compared to the 
deprojected in-plane circular velocities of all stars studied herein; these stars are mostly carbon-rich and oxygen-rich AGB stars (points). 
The red and blue points have line-of-sight velocity residuals, whose distribution is shown on the right, that exceed 50 km s _1 . These stars 
represent the tails of the distribution. 
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Figure 5. The H I gas (greyscale) and stars (points) are shown projected on the sky, where the color shading indicates the line-of-sight 
velocities after correction for the LMC's space motion. The rotation of the LMC disk is seen as a gradient running from blue in lower left 
corner to red in the upper right. ~5% of the stars, however, emphasized by a larger point size, have the opposite gradient. 
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Figure 6. Left: The color-magnitude diagram of all sources toward the LMC measured by SAGE (black points; Meixner et al. 2006). 
Points in blue are LMC stars with measured optical spectra. Green points are the stars with line-of-sight velocity residuals in excess 
of 50 km s — , which we identify as a kinematically distinct population in the LMC. Almost all of these stars appear to be AGB stars, 
with properties not grossly different from the normal LMC AGB population, with the exception that the most luminous oxygen-rich and 
carbon-rich AGB stars are not represented in the kinematically distinct population. Right: The spatial distribution of all SAGE LMC point 
sources is shown as black points. Red points show LMC AGB stars as defined by those sources that fall within the polygon in the plot on 
the left; these stars concentrate towards the LMC Bar. Blue circles show stars for which we obtained spectra, while green circles represent 
the kinematically distinct population. Although our spectra are not uniformly sampled spatially, the kinematically distinct population 
clearly avoids the Bar, contrary to the normal LMC AGB population. 
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Figure 7. The color-magnitude diagram of sources toward the SMC as measured by SAGE-SMC (black points; Gordon et al., submitted), 
plotted with four different radial cuts. The top left panel shows the SMC CMD for all sources in SAGE-SMC, while the top right, bottom 
left, and bottom right panels only include sources beyond 2°, 3°, and 4° radius, respectively, from the optical center of the SMC. The red 
points are the kinematically distinct stars seen in the LMC from this work, shifted by +0.40 magnitudes in [3.6] to account for the difference 
in distance modulus between the LMC and SMC. While the CMD containing all SMC sources has many oxygen-rich and carbon-rich AGB 
stars that are more luminous than the ones we see in the kinematically distinct population, at larger radii these luminous AGB stars 
no longer appear. There is, in particular, a good visual match between the SMC CMD at radii beyond 3° and the shifted CMD of the 
kinematically distinct population in the LMC. 



